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THE DESIGN OF A PULSE POSITION MODULATED (PPM) 
OPTICAL COMMUNICATION SYSTEM 


By Sherman Karp 
Electronics Research Center 
and 

Robert M. Gagliardi 
University of Southern California 


SUMMARY 

In recent literature the advantages of an idealized narrow 
width pulse position modulated (PPM) optical communication sys- 
tem, using coherent sources and direct photo-detection, have 
been shown. In this report, the practical design of such an 
operating PPM link is considered. System performance, in terms 
of error probabilities and information rates, are derived in 
terms of key parameters, such as power levels, number of PPM 
signals, pulse width, and bandwidths. Both background radiation 
and receiver thermal noise are included. Design procedures 
utilizing these data are outlined. Whenever possible, optimal 
design values and parameter trade-offs, in terms of maximizing 
information rate or minimizing transmitter power, are shown. 

The effect on performance of photomultipliers and their inherent 
statistics is also presented. Although the basic analysis is 
derived in terms of photon "counts," the necessary system optics 
equations are introduced to allow for overall optical hardware 
design. The primary underlying assumption is that synchroniza- 
tion is maintained at all times between transmitter and receiver. 


INTRODUCTION 

With the development of coherent sources in the optical 
region of the spectrum, there has been an increasing interest 
in the design of optical communication systems (refs. 1-6). 

The direct detection of optical radiation is presently restricted 
to photo-detection surfaces for which it has been shown that the 
released electrons obey Poisson statistics (refs. 7,8). In this 
case investigators have shown the advantages of using narrow 
width pulse position modulation (PPM) as the principal mode of 
communication (ref. 3); that is, coding information into one of 
M possible signals and transmitting it as a pulse of optical 
energy placed in one of M adjacent time intervals. It had been 
shown (refs. 1-3) that idealized versions of such systems 
optimize performance in terms of various "distance" criteria, 
and in terms of overall error probability, for cases of most 
interest. Therefore, procedures for the practical design of an 



M-ary low "duty cycle" PPM optical communication system will be 
presented in this report. In particular, performance character- 
istics in terms of key system parameters will be derived, with 
emphasis on hardware limitations and interference effects. 

Consider the PPM optical communication system shown in 
Figure 1. The transmitter is a monochromatic optical source 
operating at a fixed frequency. Information is sent by trans- 
mitting one of M signals as a pulse of optical energy at the 
same frequency, located in one of M adjacent time intervals, 
each AT seconds wide. It is assumed that complete synchroniza- 
tion between transmitter and receiver is maintained at all times; 
i.e., time coherent operation. The optical receiver detects the 
transmitted signal by attempting to determine the optical energy 
in each possible time slot, selecting the signal corresponding 
to the maximal energy. In direct photodetection this is equiva- 
lent to "counting" the number of released electrons in each AT 
interval. Background radiation entering the photodetector acts 
as erroneous energy, causing signalling errors. 


ONE OF 



BACKGROUND THERMAL 

RADIATION NOISE 


Figure 1.- Maximum likelihood processor 

In practical systems, photomultipliers are often used to 
afford an improvement in photodetection (i.e., a gain in numbers 
of released electrons) , but, unfortunately, they often behave 
randomly, which complicates system design. In addition, additive 
thermal noise may occur after photodetection. This tends to 
cause further errors in signal-decisioning . Both these latter 
effects will be considered subsequently. 

As optical radiation impinges upon a photodetecting surface, 
a series of electrons are released, each producing a current 
pulse eh(t-t m ) where e is the electron charge, t m is the time of 
release, and h(t) represents the current motion. The function 
h(t) is pulse-like, having a time width roughly equal to the 
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inverse of the photodetector bandwidth. The function h(t) is 
assumed to be identical for all electrons, and the area of h(t) 
is assumed to be normalized to unity. In the absence of thermal 
noise, the output voltage across a resistor R of the normalized 
current integrator, when sampled after AT seconds of integration, 
is then 


V 



( 1 ) 


where K is the number of electrons released during the AT second 
integration period. This result neglects "end effects," that is, 
assumes h(t) can be considered an "impulse" function with respect 
to the integrating time AT. Thus, the integrator sample is pro- 
portional to the number of electrons released in the preceding 
AT time interval, and therefore "counts" electrons. The average 
number of electrons produced in the time interval AT as a result 
of the received radiation from the transmitter will be denoted 
Kg, and is proportional to the received transmitter energy, that 
is 


K 


S 


nP R at 
hf 


(2) 


where r) is the photodetector efficiency, P R is the received peak 
optical signal power in watts, f is the mode (transmitter) fre- 
quency and h = 6.6 x 10 - 34 joules-sec. The received power can be 
related to the system optics by 


P 


R 



P t L x 10~ 4 


(3) 


where d is the optical diameter in cm, r is the range in meters, 
0 is the divergence of the antenna in radians, P T is the peak 
transmitter power, and L is the transmitting optics loss. The 
peak transmitter power can be converted to average transmitter 
power by dividing it by M, the number of signalling intervals. 

Similarly, the average number of electrons in a time 
interval AT produced by the background radiation received by the 
optical collector may be denoted by K^. Thus, 


K N = 


n p N AT 
hf 


(4) 
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where now P N is the received background radiation average power. 
This is generally written as 


P M = N, A L AX (5) 

N X R o 

where N\ is the background spectral radiance (watts/area-solid 
angle-bandwidth) , L 0 is the optical loss, Ar is the area of the 
collector, is the resolution of the receiver (solid angle) , and 
AX is the optical bandwidth. 

Note that with Pn held constant, the average number of 
"noise" electrons % is proportional to AT. This clearly in- 
dicates the advantage gained by low "duty cvcle" operation, 
i.e., using signal intervals as narrow as possible, which de- 
creases the amount of interfering radiation. The minimum value 
for AT, however, is approximately 1/AX, for then the assumption 
of fixed noise power P N is no longer valid. (For optical filters 
on the order of 5 £, AT widths of 10 - H to 10“1 2 second are 
feasible . ) 

The number of electrons counted in a signalling interval 
(i.e., an interval AT containing transmitter energy) is then a 
Poisson random variable, the average value of which is K s + %. 
For non-signalling intervals the average value is K N . Note the 
system has been tacitly assumed to be synchronized; that is, the 
integration occurs exactly during the AT data intervals. The 
maximum likelihood of processing corresponds to counting the 
number of electrons in each of the M intervals and selecting the 
interval with the largest count as the proper PPM signal. 

Allowing for likelihood draws (in which case a random selection 
is made among the drawees) , the probability of making a correct 
decision is 


P 


D 


£ 

v-=n 


1 

r + 1 


"Probability that the correct 
interval count equals r other 
interval counts and exceeds 
_the remaining M - r+1 


Temporarily neglecting the additive thermal noise and taking 
into account all the ways in which the correct interval count 
can equal r other interval counts yields: 
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Using the identity: 
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one can rewrite Eq. (6) as 
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where 


K 


a = 


N 


x-i i 
K N 


X 


* En 


i=0 


The above result is amenable to computation and can be used in 
system design to obtain performance characteristics for M-ary 
operation with fixed parameter values. An exemplary plot is 
shown in Figure 2 in which error probability P E = 1 - Pp is 
plotted versus M for fixed values of K s and K N . The results 
show the degradation in system performance as M is increased. 
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This can be attributed to the increase in likelihood draws as 
the number of intervals increases. Note that P E depends upon 
both Kg and K^, and not simply upon their ratio, so that a 
complete catalog of P E curves is required to handle all design 
conditions (ref. 9). 

Previously it was stated that the signal intervals AT 
should be as narrow as the optical bandwidths allow. This fact 
can be shown quantitatively by examining P E as a function of AT, 
when Pn is assumed to be fixed. This is shown, for example, in 
Figure 3, where M = 2 and P N is constrained such that the average 
electron noise count in an interval T Q is 10. The probability 
of error is plotted versus AT/T Q , the "duty cycle" of the trans- 
mitter. Since K s is held fixed throughout each curve, the trans- 
mitter peak power must necessarily increase proportionally, as is 
obvious from Eq. (2). Note that the error probability decreases 
monotonically as AT decreases. The minimum values at AT = 0 
are shown, but operation at these values implies infinite optical 
bandwidth. This minimum value is precisely the probability that 
a pro count occurs in signalling interval (which is also the 
noise count of the non-signalling intervals) and the receiver 
randomly selects incorrectly. It is also interesting to note 
that the pulse width AT is not particularly significant at low 
values of Kg. As Kg increases, however, the duty cycle begins 
to play a paramount role in the resulting error probability. 

Thus, attempts to increase bandwidth will have a direct payoff 
in system operation. 

It is generally tempting for communication engineers to 
base system design in terms of signal-to-noise ratios. Typically, 
this is defined as the ratio of the squared average signal elec- 
tron count to the variance of the count with background noise 
(ref. 5). In terms of a previous notation this becomes: 




K S + 


K. 


N 


( 8 ) 


To indicate the difficulty in basing design purely upon the 
signal-to-noise ratio, consider the following comparison: 

Let Kg = 10 and K N = 10; then [S/N] = 5 and, from Figure 2, 

P E = 3 x 10 - 2. On the other hand, consider the case when Kg = 5 
and K N = 0 . Again [S/N] = 5, but now P E = 3 x 10~3. Thus, with 
less signal energy, detection has been improved an order of 
magnitude with the same signal-to-noise ratio. This is due to 
the fact that Pg in Eq. (7) depends upon signal energy Kg and 
noise energy K N and not only upon the ratio. It is precisely 
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this point that distinguishes the Poisson detection problem from 
the analogous problem of transmitting one of M orthogonal signals 
over a Gaussian additive channel of equal noise power. 

It is also of interest to determine how the signal-to-noise 
ratio and error probability are related in general. By solving 
for K s in Eq. (8) one obtains: 


K, 


1 S 

2 N 


1 + 



(9) 


which can be inserted into Eq. (7). One can then plot Pg versus 
K n , for fixed values of S/N and M, as shown in Figure 4. For 
large values of K N , P^ approaches Erf (S/2N) , which is identical 
to the error probability of orthogonal coherent signalling in 
Gaussian noise with a signal-to-noise ratio of S/N. Since this 
approach is asymptotic from below, the Gaussian case will always 
be inferior to the Poisson case for the same signal-to-noise 
ratio. This point is also in agreement with results showing the 
equivalence of optimum Poisson processing and optimum Gaussian 
processing for large values of background noise power (refs. 1,4). 
Notice again that, to maintain the same signal-to-noise ratio. 

Kg must be increased. Lastly, it can be seen that even if noise 
background is negligible, K N -*■ 0 , the signalling error probability 
does not go to zero, but rather approaches 


P 


E 


V 0 



( 10 ) 


This is the same as the minimum values (at AT = 0) shown in 
Figure 3. 


EFFECT OF THERMAL NOISE AND PHOTOMULTIPLIERS 

So far in the analysis, the effect of additive thermal 
noise, which adds a random variable to the integrator output 
sample of Eq. (1), has been neglected. This complicates the 
original assertion that the receiver counts exactly the number 
of electrons at the photodetector output. Since this is the 
crux of the maximum likelihood direct detection system, it would 
be worthwhile to investigate this problem in more detail. Sup- 
pose, for example, one considers the current resulting from the 
flow of a signal photo-electron during a time interval AT = 10"^ 
second. The sample value at the integrator output across a 
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50-ohm resistor from this electron would be 8 x 10“^ volt. Now 
if the receiver operated at room temperature, the thermal addi- 
tive noise would contribute an integrator noise voltage, the 
root mean square value of which is approximately 28 x 10"® volt 
(R = 50 ohms, temp = 300°K) . Clearly, the count of a single 
electron could not easily be made in such a poor signal and noise 
condition. Photomultipliers exist, however, which effectively 
amplify the current effect of each photo-electron, resulting in 
a larger photo-electron "count" at the integrator output. Let 
this amplification factor be A, so that each photo-electron con- 
tributes a current value Ae/AT at time of sample. It would be 
interesting to determine the effect of A on P D in Eq. (7) when 
Gaussian white thermal noise of one-sided spectral level N 0 is 
added at the integrator input. If one assumes each electron 
receives the same photomultiplier gain A, the sample value due 
to the photodetector output is 


V 


KeA 

AT 


(ID 


where K is again the number of electrons produced during the 
interval AT. With Poisson statistics for the electron count 
during the AT interval, the probability density of V is then 


v v > = Z 

j = 0 


where 6 (x) is the Dirac delta function, 6 = 1/AT, and K is the 
average value of K. The thermal noise is integrated by the 
integrator and adds to the integrator sample V a random variable 
that is Gaussian-distributed with zero mean and variance N 0 B. 
Thus, the total integrator sample value z after AT seconds of 
"counting" has a probability density obtained by convolving the 
Gaussian density with the discrete density in Eq. (12), yielding: 


-K 

[jT e J 


6 (V - jeABR) 


:i2) 


P„ (z,K) 
z 


oo 



j = 0 L 


jAeBR, N q B) 


(13) 


where G(a,b,c) denotes a Gaussian density in the variable a with 
a mean b and a variance c. Observe that the sample probability 
densities are now continuous densities, and the probability of 
equal sample values occuring is zero; that is, there is a zero 
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probability of likelihood draws. Now the average count K is 
K S + K N when a signal is present in the AT interval, and is K N 
when the signal is absent. Therefore, the probability of a 
correct decision is simply the probability that the observable 
z after the correct interval exceeds the observable z after the 
M-l remaining intervals. Hence: 


z M-l 

P D = £ dZ P z (z ' K S + V [/ P z (Y 'V d ^] ‘ 


(14) 


This can be written more compactly as 


P D = 


L 

j=0 


( 


K s + K w 


- — e ^ S j dz G(z,y,c^)^ M P (z) (15) 


where 


\p (z) 


41 + 





y = AeBR 
a 2 = K 0 B 


and Erf (x) is the error integral. Equation (15) has been evalu 
ated for several values of K s and % and is shown in Figure 5, 
with 6 = 1C>9 Hz and N 0 corresponding to a noise temperature of 
300°K with a 50-ohm load. The asymptotic values for large photo 
multiplier gains are precisely the values obtained by Eq. (7) . 

At low gains, however, the thermal noise becomes the dominating 
source of error, and the probability of error increases rapidly. 
Note that to overcome the thermal environment a photomultiplier 
gain of about 10 4 is necessary for all the operating conditions 
shown. For other thermal environments one can use, as a rule of 
thumb , 

A > 600 Vtemp 


obtained by directly scaling the foregoing results. 
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In the previous analysis we have assumed that the multi- 
plier gain was a constant; that is, was the same for all photo- 
electrons. In practice, however, the gain itself is generally 
a random variable (ref. 10) with a variance, or "spread," usually 
taken as a percentage of the mean gain. It would be desirable 
now to recompute Pi-, under this situation. If one lets Aj_ be the 
electron gain of the ith released photo-electron, the integrator 
sample value for K electrons is then 

K 

V = A i e6R (16) 

i=l 


where each A^ is an independent random variable with probability 
density P A (x). The probability density of V is then obtained 
from 


P (V ) = E P (V/K) P (K) 

K=0 


(17) 


where P (V/K) is the conditional density of V, given K photo- 
electrons. Hence, from Eg. (16): 


P (V/K) 



(18) 


where ® denotes K-fold convolution. 

K 

Assume P A (x) is Gaussian with mean A and standard deviation 
aA/2 , where 0 <a < 1 represents a percentage of the gain. Then 
Eq. (18) is the Gaussian random variable with mean KAe8 and vari- 
ance K ( [aA/2]e0R) 2 so that Eq. (17) becomes: 

00 

P(V) = p- e -K G(v,iAeBR,i([aA/2]eBR) 2 ) (19) 

i=0 
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If one adds the sample contribution due to the thermal noise, 
the sample z has a probability density 


P z (z,K) = ^2 f-r e K G (v , iAegR, i ( [aA/2 ]e 3R) 2 + N 3 ) • (20) 

i=0 


Note the above equation is identical to Eq. (13), except for the 
variance terms in the Gaussian densities, and specializes to 
Eq. (13) for a = 0. Hence, the probability of detection is 
given exactly by Eq. (14), with the variance replaced by the 
above variance. The resulting error probabilities are shown in 
Figure 6 as a function of the "spreading" parameter a, using 
parameters as in Figure 5. Although the results vary somewhat 
as a function of signal and noise, one observes that for mean 
gains between 10^ and 10^ the error probabilities obtained 
earlier (Figure 2) are valid with gain spreads as high as 30 to 
40 percent. Even with spreads as high as 70 percent, the results 
indicate only about a factor of 2 increase in error probability. 
The primary conclusion then is that with suitably high mean 
photomultiplier gains, system error probabilities can basically 
be divorced of additive thermal noise effects. In this case the 
error curves plotted in Figure 2 represent the overall system 
error probabilities . 

The previous results also imply that a device average gain- 
to-spread ratio should be as large as possible for best operation. 
Consider an idealized photomultiplier characterized as a Poisson 
branching process (ref. 10). Every photoelectron emitted from 
the photoemissive surface impinges on the first stage of the 
device and releases K secondary electrons which are Poisson dis- 
tributed in number with parameter 6j_. The secondary electrons 
are then focused on the second stage where the same effect occurs. 
This process continues on through n stages resulting in a large 
electron flow at the anode for each photoelectron emitted. The 
distribution of electrons at the anode output is quite complicated 
but is unimodal and quite easily fitted with a gaussian distribu- 
tion as we have already done. In our calculations of error 
probability the two important parameters emerging were the mean 
gain of the device. A, and the mean gain to r.m.s. gain ratio, Y . 
For the idealized photomultiplier the two parameters can be 
calculated in a straightforward manner and are, in our previous 
notation , 


A 


rr « ± 


i=l 
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r 


2 


a 



6 2 6 3 


+ 


6 2 6 3 6 4 


Notice that T is almost completely characterized by the first 
stage gain 61 (a minor contribution is also made by the second 
stage) . Thus we can relate a directly to the gain of the first 
stage by 


a 



For "photon counting," Figure 6 indicates that a < .4 or that 


2 



< . 4 => 6^ > 25 . 


To take into account the effects of the remaining n-1 stages 
assume ^2 = ^3 = * * * = = ^ ' T hen 

1 = 1 - ( 1 / 6 ) n+1 < 1 = 6 
1 - ( 1 / 6 ) 1 - ( 1 / 6 ) 6 - 1 

i =0 


Therefore 6 , should be increased by 6/ (6 - 1) . Typically 6=4 
so that 



6 X > 25 • -j = 33 . 


INFORMATION RATE OF A PPM SYSTEM 

So far, only one aspect of system performance has been 
analyzed, i.e., error probabilities. The actual information 
rate that the link achieves is another important design consider- 
ation. As stated, the transmitter sends optical energy in one of 
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M time intervals AT seconds wide, thereby transmitting one of M 
possible signals in MAT seconds or at a rate log 2 M/MAT bits/sec. 
The receiver correctly determines the true signal with probability 
1 - P E and is in error with probability Pg. Because of symmetry 
the erroneous signal may be interpreted, equal-likely, as any of 
the M-l incorrect signals. Thus, the overall channel may be 
depicted as an M-ary symmetric channel in which each of the M 
possible transmitter signals is converted to itself with prob- 
ability 1 - Pg, and is converted to each other signal with prob- 
ability Pg/M-1. The information rate for such a channel is known 
to be 

_ log 2 M + P E log 2 (spf) + ( 1 _P e) l0g 2 ( 1 — P e) 

R — ~ " ' ~ v ^ -L ) 

MAT 


For convenience one may denote this as 

R = c ( K s-V M ) 

MAT 


( 22 ) 


to emphasize the dependence of the numeratpr on the stated 
parameters. By using Eq. (22) and the families of error prob- 
ability curves as shown in Figure 2, the rate R can be evaluated 
by straightforward substitution. Although specific curves for 
such a computation are not shown here, it suffices to note that 
if Kg and K N are such that P E < 10“^, Eq. (22) is to a good 
approximation , 


R 



log2 M 
MAT 


log 9 M log- M 

_ - p _ 

MAT E MAT 


(23) 


If one interprets the rate R as the source rate minus the equi- 
vocation of the channel, the PPM optical system behaves approxi- 
mately as if a source rate of log M/MAT is passed into a channel 
of equivocation Pe log M/MAT. As noted in Eq. (10), even if 
K^O (no background interference), Pg^-exp (-Kg) /2 so that the 
equivocation is not due entirely to the background noise. 

The use of Eq. (22) and the previous equations are helpful 
in determining the rate, given operating parameters. However, 
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the converse design problem, determining particular parameter 
values that achieve a desired rate, is not so straightforward. 

This is due to the fact that the rate is a somewhat complicated 
function of the parameters. Two aspects of this design problem 
that have practical application under certain operating conditions 
will be considered here: 

(1) The word period T = MAT is held fixed while the 
information bandwidth $ = 1/AT is allowed to vary 
and , 

(2) System bandwidth 3 is held fixed while the word 
period is allowed to vary. 

In both cases we are interested in the relationship between the 
rate R and the transmitter parameters Kg and M, assuming the 
noise power is held fixed. 


Fixed Work Period 

It is assumed here that AT is allowed to vary with M so as 
to maintain T = MAT constant. Thus, the system "squeezes" more 
signals into the T second period as M increases. The resulting 
rate is then 



T 


(24) 


where % T is the noise energy in T. Thus, the rate depends only 
upon the numerator of Eq. (22) . With K s fixed, increasing M 
increases the source rate, but the error probability also in- 
creases, eventiTally reaching an asymptotic value of 


-K 


P E = 


1 + K r 


L k nt " 1 + e 


NT 


-K, 


K. 


NT 


for large M. The resulting system rate increases, to within a 
constant of the entropy of the alphabet, log2 M/T . Therefore, 
it is clear that if bandwidth is expendable, one will always 
increase the system rate for large M by increasing M. in a 
practical system, this implies that one should operate with as 
wide a bandwidth as possible to exploit fully the capability of 
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the PPM system. One would, therefore, naturally consider the 
design of a system for an arbitrary rate H, when the full band- 
width (1/AT) of the system is limited. 


Fixed Bandwidth 

In this case AT is held constant (thereby fixing the noise 
energy in AT, K N ) so that both numerator and demoninator in Eq. 

(22) depend upon M, and the rate degrades quickly as M increases 
due to the log M/M dependence. A given rate, say H Q , may be 
obtained by many different combinations of Kg and M. Analytically, 
these equivalent operating points may be obtained graphically by 
noting that they are the values for which the numerator C(K S ,K N ,M), 
considered as a function of M, intersect the straight line 
H 0 ATM . By plotting these functions, for various Kg, their inter- 
section will identify (Kg,M) pairs which achieve the rate H Q . 

One may then decide on a particular operating point by invoking 
suitable design criteria. For example, one may select the 
smallest M from among the candidate pairs, which then minimizes 
the word period T = MAT. Alternatively one may choose to minimize 
the average transmitter power per information bit, which is pro- 
portional to Kg/C. In the latter case, therefore, one would 
select the operating pair (Kg,M) for which Kg/C is minimal. An 
application of this procedure is given in the next section. 


AN EXAMPLE — REAL TIME TELEVISION FROM DEEP SPACE 

To illustrate the design procedures outlined in this report, 
a television system which will transmit in real time from deep 
space will be considered. System parameter values will be chosen 
to allow use of the previously derived data. These do not always 
reflect the optimum values or current state of the art. The 
following transmission parameters will be assumed. 

Blue sky background 
Distance -- 4 x 10® km 
Optical loss — 50% 

Receiver diameter -- 16 meters (non-diffraction-limited) 
Receiver temp -- 300°K 
Optical bandwidth — 5^ 

Quantum efficiency -- 20% 

Resolution -- 1 arc sec 
Photomultiplier gain — >10~* 
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Optical frequency -- 5000 S 
Signal bandwidth — 10^ Hz. 

For real time television a rate of approximately 7 x 10^ bit/sec 
is required (corresponding to better than two samples per infor- 
mation Hz and seven bits per sample coding) . The objective here 
is to determine design parameters for a PPM system that uses 
minimal average transmitter power and a bit error probability no 
greater than 10 - 4 # Using parameters of the above table in Eqs . 
(4) and (5) yields Kjj = 0.01 as the noise background count. 
Following the discussion in the previous section, one may plot 
C(Kg,.01,M) as a function of M and determine the intersection 
with the line (H Q AT)M = (7 x 10"^)M, yielding the tabulation in 
Table I. 


TABLE I 



.8 

.9 

1.0 

1.2 

1.4 


M 

V c 

K s 

M 

V c 

19 

0.593 

1.6 

47 

0.488 

23 

0.552 

1.8 

52 

0 . 496 

28 

0.506 

2.0 

55 

0.508 

35 

0.494 

5.0 

95 

0.792 

42 

0.484 





The minimal average transmitter power occurs when Kg = 1.4, 

M = 42 which defines the PPM system. Using Eqs. (2) and (3) the 
transmitter then requires 0.072 watt. The corresponding P E can 
be obtained from Eq. (7) to get the word error probability 
(2.4 x 10~1) for which the corresponding bit error probability 
is approximately (ref. 9) P E /2 =; 10 - 1. Thus, the bit error is 
larger than that derived and, in fact, will be further increased 
by the thermal noise, as evidenced by the data of Figure 5. 

(The figure shows a slight increase in P E for Kg = 10 and M = 2 
at photomultiplier gains of 10^, and one can conclude the situa- 
tion will be worse for Kg = 1.4 and M = 42) . Thus, the minimal 
power condition is not sufficient to obtain the desired P^, for 
this example without coding. 
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To achieve the desired bit error probability Figure 2 shows 
that a system with K g = 10 and M = 100 yields bit error probabi- 
lities Pg/2 * 10“4 at the same noise level. One would expect no 
appreciable degradation from thermal noise, and Figure 6 indicates 
that only slight increases occur even with photomultipliers having 
"spreading" as high as 45 percent. The transmitter power, which 
no longer is minimal for the desired information rate, is found 
to be 0.25 watt, or an increase of 5.4 dB over the minimal average 
conditions . 


CONCLUSION 

In this paper, some design aspects of an optical M-ary 
pulse position modulated (PPM) communication system using photon 
counters (photodetectors followed by current integrators) at the 
receiver have been considered. The system considered transmits 
monochromatic optical energy in one of M time intervals, and the 
receiver determines the photon count (i.e., the received energy) 
in each interval and performs a maximum likelihood test to 
determine which signal is being received. Complete line synchro- 
nization is assumed to be maintained at all times. Performance 
characteristics are derived in terms of system parameters with 
both background radiation and thermal noise interference. In 
particular it is shown that, unlike the case of pure Gaussian 
additive noise, system performance does not depend on a few key 
parameters, but must be recomputed for different operating points. 
The important equations for deriving these characteristics are 
introduced. From these, design procedures are outlined which 
lead to best choices of transmitter power, numbers of signals, 
interval lengths, and so forth, in order to obtain desired error 
probabilities and information rates. In particular, the require- 
ments imposed upon the individual components have been indicated 
and the strong and weak points in the overall system have been de- 
lineated. No attempt has been made to estimate the cost or weight 
in building such a system. This, in fact, would be premature 
since the technology required is in its infancy and undergoing 
swift and quite radical change. In addition, no discussion of 
methods to maintain link synchronization has been included, nor 
has any consideration of atmospheric effects other than a loss 
factor been made. 
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